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Abstract. Results are presented for the 500 GeV/c pion production asymmetry and polarization of 
the Ac (Af) charmed baryon from Fermilab experiment E791. An analysis of the decay to the pKn 
final state is described. Resonant sub-channel fractions and phases are given and possible resonant 
effects in the low mass pK system discussed. Significant decay to A^ —> establishes for the 

first time the importance of a W exchange mechanism in charmed baryon decay. 


Measurements of asymmetry A,p = {d(5^ — / {d(5^ + in the yield of par- 

tiele fP and anti partiele “P ean provide information on the produetion meehanisms in¬ 
volved. Dependences of A on Xf and pp can distinguish different production models. 
Several experiments [1, and refs 1-9 therein] have shown that production of charmed 
mesons is characterized by leading particle effects and that asymmetries can be large. 
Leading particle behaviour has also been observed in production of strange hyperons in 
E791 [2] - even in a very central region. 

Branching fractions for baryon decays provide information on the relative importance 
of lowest order decay mechanisms - W exchange or spectator processes. In Ac —»■ pKn 
decay, ^ W exchange can contribute to pK*°, A*n, E*7t or pKn channels, but for the 
mode it is the only low order process possible. Evidence for this decay requires 
a large sample of pKn decays and proper analysis of interference effects in the system. 

Reported here are the first published [3] measurements of both Xf and pj dependence 
of A for charmed baryon production. We also present [4] the first full analysis of charmed 
baryon decay, measuring Ac branching fractions, relative phases and polarization. 

This study is based on a sample of 2x 10^® events produced from the interaction of 
500 GeV/c incident on thin foils, one Pt and four C. Pt target data (unequal numbers 
of n and p) were not used in the asymmetry study. The detector and data reconstruction 
are described in [5]. Cuts on geometric and kinematic quantities were made to identify 
Ac —> pK decays. The decay vertex had to be well separated (> 5a) from both 
production vertex and nearest target material. The yield, shown in Eig. l(a)-(d), was 
1,025 ±45 Ac^pK^n^ and 794±42 Ac—^pK^n^. Events were divided into 5 regions 
of Xf and 5 of in the overall ranges: —0.1 <Xf < 0.6 and < 8 (GeV/c)^ chosen 
to have clear Ac signals in each. Eits similar to those shown in the figure were made to 
each sample to determine the number N{Ac) and A^(Ac) of signal events in each range. 


^ Representing the E791 Collaboration 

^ Note that charged conjugate states are implied unless stated otherwise. 




FIGURE 1. Ac(Ac) samples: (a) A^- {xp < 0); (b) A^- {xp < 0); (c) A^ {xp > 0); (d) Ac {xp > 0). Fits 
shown in (a)-(d) are Gaussian peaks on 2nd order polynomial backgrounds. Asymmetries: (e) vs xp (f) vs. 
Pj. The solid curves in (e) and (f) are the prediction of Pythia/Jetset. 


Ac/Ac ASYMMETRIES 

Efficiencies £ (e) for Ac (Ac) were not quite equal due to the asymmetric effect of the 
intense beam on the drift chambers. This effect was greatest at large xp and low p^. 
It was necessary therefore to estimate the ratio r = e/e in each of the 5 xp and 5 pj 
ranges using Monte Carlo samples of Ac & Ac generated with Pythia/Jetset, projected 
through a simulated E791 detector and subjected to the same reconstruction code and 
selection criteria as the data. Corrected asymmetries A = (N — N/r)/{N + N/r) were 
then obtained in each range. 

The main sources of systematic uncertainty (parametrization of signal and back¬ 
ground shapes and precision of r) amounted to less than 50% of the statistical uncer¬ 
tainty in all instances. The results are shown in Eigure 1(e) and (f) and compared with 
earlier studies in Table 1. 

TABLE 1. Comparison with asymmetries (%) from earlier k^N experiments. 
xp region E791 ACCMOR [6] SELEX [7] 

Xf<0 20±10±6 

xp>0 12.3±3.7±1.6 0.5±7.9 25±15 


The asymmetry is positive and flat throughout the range. This might result from 
the additional energy required to produce additional baryons when a Ac is produced, 
favouring Ac production in general. The solid curve in Eigure 1 is the prediction of 
Pythia/Jetset and clearly does not describe the data well. Two component intrinsic 
charm/coalescence models [8], [9] predict a rising asymmetry beginning at the low end 
or possibly below the range of this data. Eeading particle effects would also result in a 
























rising asymmetry in the entire Xf <0 region. The data do not rule that possibility out. 


ANALYSIS OF THE DECAY Ac ^ PK-n+ 


A eleaner sample was required for this analysis. The length cut was increased to 8a 
and a neural net criterion was used to optimize the significance S / ^/S + B of the signal 
(S = 886 ± 43) over background (B ~ 300) in the fit region. 

These decays were defined by five independent variables, e.g. two Dalitz plot coor¬ 
dinates and orientation of the decay plane relative to production plane (z axis). The Ac 
could have polarization 

Each isobar decay channel Ac~^R{^ab)c was assigned an amplitude labeled by the 
z component of Ac spin, m, and proton helicity (±^) in the Ac rest frame 
A* i_ = 5*(Mflfe)(a±e'“±|m, >) where are 

in,±2 

the two possible helicities for R with unknown coefficients & b±e'^^. The 

B’^{Mab) were Breit Wigner functions. For non-resonant NR decay to pKn a similar 
amplitude with B = \ was used. An unbinned, maximum likelihood fit was used to 
determine a±, a±, Z7±, (3± and three values for (one in each of three xp ranges). The 
signal probability density function was 


Ts = ^X£X[(l-fEA 
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Five dimensional efficiency (£) and background density were estimated empirically from 
a MC sample and Mpp-n sidebands. Modes included were pA'*°(890), A++(1232)^^, 
A(1520)7t+ and AT?. 

The fit shown in Figure 2 is seen to be good except for the low mass K p region where 
an unmodelled enhancement is seen. Many T* exist which could possibly account for 
this. Each Y*k channel added to our fit requires > 4 more parameters making it difficult, 
with our limited sample, to include more than one Y*. Adding A(1600)7t, E(1600)7t or 
the tail of the E(1405)7t alone made no significant improvement. 

Isobar fractions fp were computed by integrating over the five dimensions of the fit x: 
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Branching ratios with respect to pK^n^ are compared with earlier results in Table 2. 

Good agreement is seen, but the significance of signals from NA32, where only mass 
projections were fit is overestimated. NA32 errors are comparable to E791, but E79rs 
sample is much larger. This is because correlations among channels and relative phases 
were neglected in the NA32 analysis. 

The mode is comparable to pK and clearly significant. 
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FIGURE 2. Left: Projections of fit (dashed lines) onto three mass pair and three angular variables (solid 
lines). Data lie in the range 2265 < M{pK^ 7t+)< 2315 MeV/c^. Right: A^/Ae Polarization from fit. 


TABLE 2. First three columns are branching fractions relative to total pKn mode (% corrected for 
unseen decays). The last four columns are resonant phases, described in the text, measured only by E791. 

Mode E791 NA32 ISR E791 relative phases (degrees) 


[10] [11] a+ p+ a 


pK*° {S90) 

29±4±3 

35+® +3 

42 + 24 

58 + 28 

135 + 38 

198 + 24 

303 + 32 

A++(1232)A- 

18 + 3 + 3 

123 ±5 

40+17 

285 + 23 

280 + 23 

= a+ 

= P+ 

A(1520)7t+ 

15+4 + 2 

9+^ + 2 

- 

340 + 30 

-3 + 32 

= a+ 

= P+ 

NR 

55 + 6 + 4 

561^ + 5 

- 

199 + 31 

0 (fixed) 

43+41 

65 + 21 


SUMMARY 

A^. production asymmetry in the range —0.1 < xp < 0.6 and < 8 [GeY jc)^ is 
eonstant at ~ +0.15 favouring Ac over Ac. Models requiring a rising asymmetry toward 
negative xp are not ruled out however. An amplitude analysis of the Ac deeay shows the 
Ac —> mode to be large indieating that the VL exehange amplitude is important. 
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